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Abstract — Energy harvesting is a promising technique to 
overcome the limit on energy availability and increase the 
lifespan of battery-powered embedded systems. In this paper, the 
question of how one can achieve the prolonged lifespan 1 of a 
real-time embedded system with energy harvesting capability 
(RTES-EH) is investigated. The RTES-EH comprises a 
photovoltaic (PV) panel for energy harvesting, a supercapacitor 
for energy storage, and a real-time sensor node as the embedded 
load device. A global controller performs simultaneous optimal 
operating point tracking for the PV panel, state-of-charge (SoC) 
management for the supercapacitor, and 

energy-harvesting-aware real-time task scheduling with dynamic 
voltage and frequency scaling (DVFS) for the sensor node, while 
employing a precise solar irradiance prediction method. The 
controller employs a cascaded feedback control structure, where 
an outer supervisory control loop performs real-time task 
scheduling with DVFS in the sensor node while maintaining the 
optimal supercapacitor SoC for improved system availability, 
and an inner control loop tracks the optimal operating point of 
the PV panel on the fly. Experimental results show that the 
proposed global controller lowers the task instance drop rate by 
up to 63% compared with the baseline controller within the same 
service time (i.e., from sunrise to sunset.) 

Index Terms — Energy harvesting, real-time embedded system, 
photovoltaic panel, supercapacitor. 



I. Introduction 

E NERGY consumption is one of the most important design 
criteria for battery-powered devices and systems. 
Dynamic power management (DPM) [2] [3] [4] [5] and dynamic 
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voltage and frequency scaling (DVFS) [6][7][8][9] techniques 
have been extensively investigated for reducing energy 
consumption of battery-powered systems while meeting 
performance targets. However, even with sophisticated DPM 
or DVFS techniques, batteries still need to be charged or 
replaced for sustained system operation. 

There are some applications that do not allow battery 
charging or replacement, such as sensor nodes deployed in 
radioactive surroundings. Therefore, energy harvesting 
techniques have been deployed to overcome the limit of 
energy source and increase the lifespan of untethered 
embedded systems [10][1 1][12]. Energy harvesting from 
ambient energy sources including sunlight, wind, and tidal 
wave can provide unlimited and inexhaustible power supply 
for embedded systems. However, the harvested energy 
exhibits intermittency characteristics. For instance, the output 
power of a photovoltaic (PV) panel at noon is an order of 
magnitude higher than that in the morning or evening. 
Therefore, an energy storage element, typically a rechargeable 
battery [10] or a supercapacitor [12], is mandatory in an 
embedded system with energy harvesting capability for 
sustainable power supply. Several prototypes have been 
presented such as Heliomote [10] and Prometheus [11]. An 
interesting summary of various design considerations for an 
embedded system with energy harvesting is provided in [13]. 

A number of recent work have focused on power 
minimization algorithms and techniques for a real-time 
embedded system with energy harvesting capability 
(RTES-EH) [14] [15] [16]. These work are different from the 
traditional real-time task scheduling and DVFS work 
[6] [7] [8] [9] because of the intermittent nature of harvested 
energy and limited size of the energy storage element. Indeed, 
the primary purpose of the RTES-EH task scheduling is to 
enhance the system availability in spite of the sporadic nature 
of the harvested energy. The lazy scheduling algorithm 
proposed in [14] executes tasks as late as possible but at full 
speed. DVFS combined with lazy scheduling is used [15] [16] 
to reduce deadline miss rate minimizing energy dissipation. 

In this work, we present a global control algorithm for a 
RTES-EH comprised of a PV panel as the energy harvesting 
source, a supercapacitor as the energy storage element, and a 
real-time sensor node as the embedded load device. The 
presented global controller performs simultaneous optimal 
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operating point tracking of the PV panel, state-of-charge 
(SoC) management of the supercapacitor, as well as real-time 
task scheduling with DVFS for the sensor node. The combined 
task scheduler and DVFS manager in turn relies on a precise 
solar irradiance prediction method. 

Limitations of the previous work and key contributions of 
the present work are summarized in the following. 

1) Using battery as the energy storage element in previous 
work [14] [15] [16] may prevent a RTES-EH from having 
a long lifespan. Typically, batteries have a rather short 
cycle life (up to several hundred charge-discharge cycles) 
after which they have to be replaced. Alternatively, 
supercapacitors constitute a highly effective energy 
storage element for a RTES-EH due to (i) 
orders-of-magnitude longer cycle life than batteries, (ii) 
higher power capacity to deal with bursty power demands 
of load devices, and (iii) very low environmental impact 
[17]. Because the terminal voltage variation of a 
supercapacitor is much larger than that of a battery [18], it 
is not possible to apply previous work directly to the 
supercapacitor-based RTES-EH. In this work, we present 
a global control algorithm for a supercapacitor-based 
RTES-EH, accounting for supercapacitor characteristics. 

2) The converters connecting the energy source to the energy 
storage element, and the energy storage element to the 
load device cause power loss in a RTES-EH. The power 
loss is significant in a supercapacitor-based RTES-EH due 
to large variation of the supercapacitor terminal voltage. 
The present work considers the converter efficiency 
variation based on accurate power converter modeling in 
order to optimize the overall system availability (this is in 
clear contrast to previous work which assumes 100% 
power conversion efficiency). 

3) Most of the previous work assumes ideal component 
models for both the energy harvesting source and the 
energy storage element. The present work employs 
accurate component models to take into account the 
nonlinear current-voltage (I-V) characteristics of the PV 
panel, and energy loss and self-discharge of the 
supercapacitor when optimizing the system availability. 

4) To derive a global control algorithm for the RTES-EH 
with a PV panel as the energy harvesting source, it is 
necessary to predict the solar irradiance profile, and 
thereby, estimate the PV panel energy generation profile. 
The present work uses an accurate solar irradiance 
prediction algorithm, which refines the solar irradiance 
prediction at each decision epoch based on actual solar 
irradiance values observed previously. 

To develop an effective global control algorithm, there are 
two fundamental subproblems to be solved, i.e., the 
supercapacitor charging problem and the converter-aware 
frequency setting problem. The first problem aims at 
dynamically maximizing the supercapacitor charging current 
by optimally setting the operating point of the PV panel (this 
is similar to the problem discussed in [18] [19]). Different from 
the traditional maximum power point tracking (MPPT) 
techniques [20], which only focus on maximizing the output 



power of the PV panel itself, the supercapacitor charging 
problem accounts for the power loss in the converter to 
maximize the input power to the supercapacitor. The second 
problem focuses on finding the optimal sensor node execution 
frequency for a single task instance without deadline 
constraint such that the amount of energy extracted from the 
supercapacitor is minimized. According to [21], the minimum 
execution frequency in the sensor node is no longer always 
optimal in energy saving for a task instance without deadline 
constraint when the converter power loss is taken into account. 

Based on the optimal solutions to the aforesaid two 
subproblems and the PV panel energy generation profile, we 
present a near-optimal global control algorithm by using a 
cascaded feedback control structure. An outer supervisory 
control loop maintains the supercapacitor SoC throughout the 
system operational period (from sunrise to sunset) so as to 
enhance the overall system availability and reduce the task 
instance drop rate. The supercapacitor SoC maintenance is 
achieved through effective real-time task scheduling combined 
with DVFS as well as selective task instance dropping in the 
sensor node. An inner control loop, on the other hand, tracks 
the optimal operating point of the PV panel (and hence the 
input current of the supercapacitor) on the fly to maximize the 
amount of harvested energy. Experimental results demonstrate 
that the proposed accurate component model-based global 
control algorithm on the RTES-EH significantly lowers the 
task instance drop rate by up to 63% compared with baseline 
control algorithm. 

The remainder of this paper is organized as follows. Section 
II presents models of components in the RTES-EH. Section III 
introduces the RTES-EH architecture and the real-time task 
set. Section IV defines a formal problem formulation of the 
RTES-EH global control problem. And Section V provides the 
optimal solutions to the aforesaid subproblems. Section VI 
discusses about the solar irradiance prediction algorithm. 
Section VII proceeds with the global control algorithm. 
Section VIII and IX are experimental results and conclusion, 
respectively. 

II. Component Models 

A. PV Cell Model and Characterization 

As stated before, the energy harvesting source in the 
RTES-EH is a PV panel. The PV panel consists of a number 
of PV cells connected together in a balanced series/parallel 
configuration i.e., an mxn array, where m is the number of 
PV cells connected in series and n is the number of PV cells 
connected in parallel. We use V pv and I pv to denote the output 
voltage and current of the PV panel, respectively, and V pvc 
and I pv c to denote the output voltage and current of a single 
PV cell inside the panel, respectively. The output voltages and 
currents have the following relationships: 

V pv — 771 ' V pV)C , l pv — 71 ' I pv c . (1) 
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Fig. 1. An equivalent circuit model of a PV cell. 



by the micro-controller. Fig. 2 shows a conceptual structure of 
the PWM buck-boost power converter. We denote the input 
voltage, input current, output voltage, and output current of 
the power converter by V in , I in , V out , and I out , respectively. 
The power loss of the converter, P conv , consists of the 
conduction loss, the switching loss and the controller loss [23]. 

According to the energy conservation law, P con v satisfies: 



A typical equivalent circuit model of a PV cell is shown in 
Fig. 1, with I-V characteristics given by: 

Ipv,c — h — Id — Ish (2) 

= 7 l (G) - I 0 (T) ( e (V.c+W R s) ' j h - l) - V- c+ W-^ 

V / Rp 

where 



F in ' I in Pconv F F out ' I out- (^) 

Based on the relationship between V in and V out , the power 
converter operates in one of the two possible operating modes: 
the buck mode when V in > V out and the boost mode, 
otherwise. When the converter is operating in the buck mode, 
its power loss P conv is given by 



and 



k(G) — - Il(G ST c), 
Gstc 



Pconv - tout (^L "F H R swl + (1 - H) • R sw 2 + R sw 4 ) (6) 

(3) (A/) 2 

+ --(R L + H -R swl + (1-H)- Psw 2 F Rsw4 + Pc) 

F Fin ’ fs ' ( Qswl F Qsw2) F Fi n • I controller? 



3 QEg f 1 _ l\ 

Io(T) = IoVstc) ■ (^) ■ e ^ stc — F/ e (4) 

Parameters in (2) - (4) are defined as follows. G denotes the 
solar irradiance level; T is the cell temperature; q is the charge 
of the electron; E g is the energy bandgap; and k is the 
Boltzmann’s constant. STC stands for standard test condition 
in which the irradiance level is 1000 W/m 2 and the cell 
temperature is 25 °C. There are still five unknown parameters, 
which are typically not provided by manufacturers, yet to be 
determined. These unknown parameters are: 

• Il(G S tc) : photo-generated current at STC, 

• I 0 (T stc ): dark saturation current at STC, 

• R s \ PV cell series resistance, 

• R p : PV cell parallel (shunt) resistance, and 

• A: ideality factor of the diode. 

We adopt the method proposed in [22] to extract the 
above-mentioned five unknown parameters from the measured 
PV cell I-V curve at any specific environmental condition. 
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Fig. 2. PWM buck-boost power converter architecture. 



B. Power Converter Model 

Both the charger and the DC-DC converter used in the 
RTES-EH are programmable PWM (pulse width modulation) 
buck-boost power converters. A charger regulates its output 
current to the value set by the micro-controller, while a 
DC-DC converter regulates its output voltage to the value set 



where H = V out /V in is the PWM duty ratio and A / = V out • 
(1 — H)/(Lf • f s ) is the maximum current ripple; f s is the 
switching frequency; I controller * s the current flowing into the 
micro-controller of the power converter; R L and R c are the 
internal series resistances of the inductor L and the capacitor 
C, respectively; R swi and Q swi are the turn-on resistance and 
gate charge of the z-th MOSFET switch shown in Fig. 2, 
respectively. 

On the other hand, the converter power loss P con v h 1 the 
boost mode is given by 

Pconv = fe) 2 • (7) 

(Rl + H ' Psw3 F (1 H)R sw4 + Pswl + H( 1 - H)R C ) 

+ ^f(R L + HR SW 3 + (1 - HXRsw* + Rc) + Rswi) 

FF 0Ut • f s • (Q sw 3 + Qsw4) F V in • I controller ? 

where H = 1 - V in /V out and A / = V in • H/(L f • f s ). 

It can be observed from (6) and (7) that the converter power 
loss P CO nv i s a function of V in , F out , and l out . Notice that P conv 
can also be written as a function of V in , I in , V out , which is not 
provided in this paper due to space limitation. 



C. Supercapacitor Model 

Supercapacitors have a very small internal resistance and 
thus very small IR power loss is incurred during charging and 
discharging. Supercapacitors also exhibit a significantly higher 
volumetric power density and a longer cycle life compared 
with batteries [17]. These features make supercapacitors a 
good choice for energy storage element in a RTES-EH. 

A primary disadvantage of supercapacitors is the full-range 
terminal voltage variation as a function of its SoC. The 
terminal voltage of a supercapacitor, denoted by V cap , is a 
linear function of its SoC denoted by Q cap • More precisely, 
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Kap = Qcap/Ccap = pE cap /C cap , where E cap is the energy 
stored in the supercapacitor and C cap is the capacitance of the 
supercapacitor. The terminal voltage variation of a 
supercapacitor leads to power conversion efficiency variation 
in the RTES-EH, which motivates us to address this issue in 
the proposed global control algorithm. 

Moreover, a supercapacitor may lose as much as 40% of its 
energy per day due to self-discharge [17]. The terminal 
voltage decay of a supercapacitor after a time period At is 
given by: 

_At 

V cav (t + At) = V cap (t) ■ e r, ( 8 ) 



denote the output voltage and current of the PV panel at time t 
by V vv (t) and/ pi 7 (t), respectively. We control the PV panel 
operating point (V pv ( t),I pv (t)) through effectively setting the 
charger output current l C ap,in(P) as shown in Fig. 3. 

We denote the terminal voltage, input current and output 
current of the supercapacitor at time t by V cap {f), I cap j n (t), 
and I cap, out (.ty 5 respectively. The supercapacitor energy 
E ca V (t) at time t is given by 

E cap O') — E cap (f sunrise) Jr ^ capO ) \^cap,inO) ~ 

^cap.outO) ~ Icap.sd (t)) dr, (9) 



when no external power supply or load is connected to the 
supercapacitor. The parameter T in ( 8 ) denotes the 
self- discharge time constant. 

III. Real-Time Embedded System with Energy 
Harvesting Capability 

A. System Architecture 

This paper investigates a RTES-EH with the architecture 
shown in Fig. 3. In particular, the RTES-EH is comprised of a 
PV panel as the energy harvesting source, a supercapacitor as 
the energy storage element, a real-time DVFS-enabled sensor 
node as the load device, and power converters connected in 
between. A charger connects the PV panel to the 
supercapacitor, and a DC-DC converter connects the 
supercapacitor to the real-time sensor node. A 
micro-controller controls the system operation and executes 
the control algorithm of the system. 



^Micro-Controllerj 



Energy Harvesting Module Embedded Load Device 




Fig. 3. Block diagram of the RTES-EH. 



The RTES-EH operates from sunrise T sunrise until sunset 
T s Un set the system operational period) and it is in idle 

state for the rest of time during a day. We denote the solar 
irradiance and temperature on the PV panel at time t by G (t) 
and 7(t) , respectively. Since the temperature has only 
secondary effect on the I-V characteristics of the PV panel, we 
do not consider the PV panel temperature change in this paper, 
i.e., we assume constant PV panel temperature 7(t) = Tstc • 
On the other hand, G(f) has a huge impact on the I-V 
characteristics of the PV panel. We thus estimate the G(f) 
profile by using a novel and highly accurate solar irradiance 
prediction method as shall be discussed in Section VI. We 



where E C ap(T sunrise ) is the initial supercapacitor energy at 
time T sunrise , and/ capsd (t) is the self-discharge current at 
time t, which is calculated based on (8). 

We assume a single-core microprocessor as the sensor node 
in the RTES-EH. We denote the supply voltage, the input 
current, and the power consumption of the sensor node by 
Vpr OC (t)i IprocO') andP proc (t), respectively. The sensor node 
has N discrete operating frequencies: freq min = freq 1 < 
freq 2 < ••• < freq N = freq max ; each operating frequency 
freq n corresponds to a supply voltage V proc>n , an input current 
Iproc.m an d a power consumption value P proc>n (1 < n < N). 

The input voltage, input current, output voltage, and output 
current of the charger are V pv (f) , I pv (f) , V cap (t ) , and 
I cap,in(t)i respectively. The charger power loss at time t, 
denoted by P CO nv ,inO)> is a function of V pv (t), I pv (t), and 
V cap (t) as discussed in Section II.B. Similarly, the power loss 
of the DC-DC converter, denoted by P CO nv,out0 ), i s a function 
of its input voltage, output voltage, and output current, i.e., 
K:ap(^)? Vproc(P) ? an d / proc (t), respectively. The following 
two equations hold due to the energy conservation law: 



V pv (t)i pv (t) = 

Vcap(.t)Icap,in(t) Pconv, in (t), (io) 
Vcap(t)Icap,out(t) — f p roc(^')^proc(^') Pconv, out ) • (H) 



B. Real-Time Task Set 

In this paper, we assume a frame-based, preemptive, hard 
real-time task set in the sensor node [ 8 ]. The set of M 
real-time periodic tasks is denoted by [T lt T 2 , T M }- Task 
instances of T t (1 < i < M) are released periodically with a 
period of D t , which is also the relative deadline of task 
instances off. Task instances of f have a workload ofVl^, 
which is given as the number of clock cycles required to 
complete a task instance. We use the independent-task 
assumption that tasks (and therefore task instances) are 
independent from each other [ 1 3] [ 14] [ 1 5] . 

We define the frame length as Plcm j which is the 
hyper-period of all the tasks (i.e., the least common multiple 

of D lt D 2 ,-" ,D m ). We use K = Tsunset T sunrise tQ d eno t e the 

t lcm 

total number of frames from T sunrise to T sunset . Therefore, 
each k- th frame (1 < k < K) begins at time T sunrise + (k — 
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1) ■ T lcm and ends at time T sunrise + k ■ T LCM . Fig. 4 shows 
the start time and end time of the k- th frame and also the time 
windows of task instances in the frame. The duration of a 
frame is in the order of minutes or seconds, which is much 
shorter than the RTES-EH operational period. We use 5 
minutes for Tlcm in the simulation to represent a reasonable 
trade-off between computational complexity and control 
accuracy of the global control algorithm. 



( y rj->2,k 

Task 1 |* — - — ^ — *- 



A:-th frame 




1 T 

Task i 'f — L 



i,k 



T 2,k 



rn\,k 



rr>2,k 



Task M 'f- 

T sunrise + _ 1) ' ^ LCM 

Fig. 4. Task instances in the k- th frame. 
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The earliest deadline first (EDF) scheduling [6] is 
incorporated in the sensor node to schedule the periodic tasks. 
The y-th (1 < j < T LCM /D i ) task instance of T t in the k- th 
frame, denoted by T^ ,k , is released at time T sunrise + 
(/c — 1 )T lcm + (J — 1 )Dj, which is the absolute release time 
of that task instance. Hence, the absolute deadline of task 
instance T l ’ ,k is equal to T sunrlse + (k - 1 )T LCM + jD t . 



IV. Problem Formulation 

We use fsfit) to denote the sensor node operating 
frequency assigned for executing an instance of T t at time t. 
For real implementation, we assume that fsfit) can only take 
discrete values in the set {freq lf freq 2 ,--- , freq N }. If the 
sensor node does not execute any task instances of T t at time t, 
we have fsft ) = 0 . The single-core microprocessor 
assumption of the sensor node ensures that at most one of the 
fsft) values for 1 < i < M at time t is non-zero. Therefore 
the operating frequency of the sensor node 
fs(t) = Hi<i<M/ 5 i(0- We adopt a general definition of a 
schedule of task instances where task instance dropping or 
timing constraint violation is allowed to accommodate the 
intermittent nature of the harvested energy of the RTES-EH 
[15][16]. We define a new variable ST k such that sT k = 1 
indicates the completion of task instance Tf before its 
deadline, i.e., 



fT sunrise +(.k-l)T LC M+j D i 
^sunrise + (k — U^LCM + (/ — 1 i 



fsi(t)dt > W t . 



( 12 ) 



On the other hand, ST k = 0 indicates the task instance 
has missed the deadline, i.e., 



fT sunrise +(.k-l)T LC M+j D i 
Tsunrise + (k ~DTlcM + 0” — 4)U j 



fsi(t)dt < W t . 



(13) 



We minimize the task instance drop rate, or equivalently, we 
maximize total number of fully executed task instances, given 
by 



T otal_Executed — Yii<k<K Hi<£<m ^i<j<T LCM /Di $1’ ■ (14) 

The power consumption P proc (t ) and supply voltage 
V p roc (0 °f the sensor node are dependent on fs(t). The 
supercapacitor discharging current I ca p,out (0 I s determined by 
Pproc(f) , P C onv,out(t ) » and V cap (0 > using (11). The 
supercapacitor charging current l C ap,in(t) I s determined by 
(Vpv(t)>lp v (t)), P CO nv,m(t)> and v cap (t), using (10). We 
calculate E cap (t ), i.e., the energy stored in the supercapacitor 
at time t, using (9). Therefore, the fsft) profiles for 
1 < i < M and the PV panel operating point (V pv (t), I pv (f)) 
are the control variables of the RTES-EH. We achieve the 
supercapacitor SoC control (i.e., the control of its energy) 
through task scheduling with DVFS (i.e., finding th q fsft) 
profiles) and PV panel operating point tracking. 

The formal statement of the global control problem for the 
RTES-EH is given as follows: 

RTES-EH global control problem statement: 

Given the initial supercapacitor energy E cap (T sunrise ), 

Find a schedule and frequency allocation of all tasks, 
represented by fsft) ( 1 < i < M ), and the PV panel 
operating point (V pv (t), I pv (t )) for t E [T sunrise 

> Tsunset\’ 

Maximize the total number of fully executed task instances by 
the deadline, as given by (14), 

Subject to the energy conservation law (9), (10), and (11). 

Please note that the solar irradiance profile G(t) is not 
given explicitly in the problem statement. Instead, G(t) will 
be estimated by a prediction method in the global control 
algorithm. In the next three sections, we will introduce two 
optimization subproblems and their solutions, the solar 
irradiance prediction algorithm, and the global control 
algorithm based on the solutions to the subproblems and solar 
irradiance prediction. 



V. Optimization Subproblems 

We first define two fundamental optimization subproblems 
that are necessary in solving the original global control 
problem. The two subproblems are the super capacitor 
charging (SC) problem and the converter- aw are frequency 
setting (CA-FS) problem. These two subproblems are inherent 
from the architecture of the RTES-EH as shown in Fig. 3. 
According to the energy flow in the RTES-EH, the RTES-EH 
can be partitioned into two parts, i.e., the energy harvesting 
part and the energy consumption part. The energy harvesting 
part consists of the PV panel, the supercapacitor and the 
charger in between. The energy consumption part consists of 
the supercapacitor, the sensor node and the converter in 
between. For the energy harvesting part, one natural question 
is how to collect the maximum amount of energy from the PV 
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panel to the supercapacitor, and the SC problem answers the 
question. For the energy consumption part, one natural 
question is what is the optimal sensor node frequency such 
that the energy extracted from the supercapacitor is minimized, 
and the CA-FS problem answers the question. 

A. Supercapacitor Charging (SC) Problem 
The SC problem aims at maximizing the supercapacitor 
charging current 7 cap j n (t) by optimally setting the operating 
point (V pv (f), I pv (t)) of the PV panel, given the current solar 
irradiance G(t ) and the amount of energy stored in 
supercapacitor E cap (t) . Solving the SC problem optimally at 
any time t E [T sunrise> T sunset \ is required for the optimal 
solution to the global control problem for the RTES-EH. The 
SC problem is nontrival and requires close examination of 
characteristics of both the PV panel and the charger. 




Fig. 5. I-V and P-V characteristics of the PV panel under different solar 
irradiance values. 



Fig. 5 illustrates the I-V and power- voltage (P-V) 
characteristics of the PV panel under different solar irradiance 
G values. For a specific G value, the PV panel output current 
I pv decreases as the output voltage V pv increases, and the PV 
panel output power V pv ■ I pv is maximized when it operates at 
the maximum power point (MPP), marked by red dots in Fig. 
5. Conventional MPP tracking (MPPT) technique maintains 
the PV panel operating point at its MPP. However, such 
technique cannot guarantee the maximum amount of energy 
transferred into the supercapacitor due to the neglect of 
efficiency variation of the charger [18]. Recent maximum 
power transfer tracking (MPTT) technique takes into account 
the non-ideal behavior of the charger, thereby maximizing the 
energy transferred into the supercapacitor [18]. In this paper, 
we adopt the MPTT technique to solve the SC problem 
optimally at any time t € [T sunrlse , T sunset ], 

There are two possible implementations of the SC problem 
solution. The first one is software-based or model-based 
implementation whereas the second one is hardware-based 
implementation. The former implementation uses accurate PV 
panel modeling as described in Section II.A and derives the 
optimal supercapacitor charging current 7 cap j n (t) . This 
implementation is adopted in the supervisory control of the 
RTES-EH. On the other hand, the latter implementation is 



based on a perturb and observe (P&O) heuristic in a hardware 
micro-controller to dynamically track the optimal operating 
point of the PV panel. This implementation is adopted in the 
inner control loop of the RTES-EH and does not need prior 
knowledge of the PV panel modeling. 

B. Converter-Aware Frequency Setting (CA-FS) Problem 
The CA-FS problem aims at finding the optimal sensor 
node frequency for executing a specific task instance such that 
the energy extracted from the supercapacitor is minimized. 
Recall that executing a task instance at minimum clock 
frequency is not necessarily the minimum energy solution [21], 
because the power conversion efficiency of the DC-DC 
converter in the system is a function of the terminal voltage of 
the supercapacitor, which in turn is dependent on the amount 
of energy stored in the supercapacitor E cap (f). 

We derive the total energy extracted from the 
supercapacitor during the execution of a task instance as a 
function of the execution frequency /. Note that / can only 
take discrete values from the set [freq lf freq 2 w , freq N }. 
The supply voltage and power consumption of the sensor node 
during execution of a task instance are functions of/, i.e., 
Vproc(f) and P proc (/), respectively. The total execution time 
of a task instance off, denoted by T t exe (f), is given as 
T? xe (f) = W t /f, where W t is the workload of task instances 
off as discussed in Section III.B. We assume that the V cap 
value will not change significantly during the execution of a 
single task instance, since the execution time is too short to 
make noticeable SoC change in the supercapacitor. The power 
loss of the DC-DC converter during execution of a task 
instance is a function of its input voltage l^ ap , output voltage 
VprodD, and output current I proc (f) = P pr oc(f)/ V proc(f)- as 
specified in Section II. B. Therefore, we derive the power loss 
of the DC-DC converter as a function of/andF cap , i.e., 
Pconv.outif >Vcap)- The amount of energy extracted from the 
supercapacitor for executing a task instance off, denoted by 
Extract ( f>Vcav) isgivenaS 

E fxtract (fiVcap) = T exe (f) . + P conv>out (f ,V cap )), 

(15) 

which is not a monotonically increasing function of /. This 
implies that using the minimum clock frequency for task 
execution does not always minimize the total energy drawn 
from the supercapacitor. We denote the optimal sensor node 
frequency for executing a task instance of f by f° vt (y cap ) 
that minimizes the total energy drawn from the supercapacitor. 
fi° Pt (Vcap) is calculated as 

fi V (Vcap) = ^ r & m ^ n freq min <f<freq max Ei (/ >^cap) • (16) 

VI. Solar Irradiance Prediction Algorithm 
The solution to the global control problem relies on the 
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solar irradiance prediction results. In this section, we introduce 
a solar irradiance prediction algorithm based on [24]. From the 
start time T sunrise to the end time T sunset of the RTES-EH 
operational period, there are a total number of K frames of 
task instances to be executed, as mentioned in Section III.B. 
We assume that the solar irradiance is constant during a frame. 
This assumption is valid since the duration of a frame is in the 
order of minutes or seconds. Therefore, we will predict the 
solar irradiance for rest of the day at the start time of the k- th 
(for 1 < k < K ) frame in the d- th day, i.e., T* k , where 
subscript s denotes “start”. 

For the solar irradiance prediction, an important observation 
is that the actual solar irradiance over the k- th frame in the 
d- th day may be viewed as the solar irradiance over the k- th 
frame in a sunny day multiplied by a decay factor , 
representing the effect of clouds, if the d- th day is cloudy. 
Obviously, the sunny day solar irradiance over the k- th frame 
varies with season changes. This effect is however captured by 
a smoothing operation as shall be discussed later. In general, 
the solar irradiance prediction algorithm consists of the initial 
prediction and the intra-day refinement: the initial prediction 
is performed at the beginning of each day i.e., T* fl to predict 
the sunny day solar irradiance over each frame during that day; 
the intra-day refinement is performed at the start time of the 
k- th (for 1 < k < K ) frame in that day to predict the decay 
factors and thereby refining the prediction results over the rest 
of the frames in that day. Please note that at the start time of 
the k- th frame, the decay factor in the (/c — l)-st frame is 
already known. And decay factors of different frames in the 
same day are positively correlated. 

In the initial prediction phase, we adopt a variant of the 
well-known exponential average-based prediction method, in 
order for effectively predicting the sunny day solar irradiance 
in each frame of a day (please see the discussion below for the 
differences between our version and the standard method.) 
Consider that we are at time T* ± = T sunrise of the d- th day. 
We want to derive the prediction value of the sunny day solar 
irradiance in the k- th (for 1 < k < K ) frame of that day, 
denoted by Predict k , based on the prediction value of the 
sunny day solar irradiance in the k- th frame of the (d — l)-st 
day, denoted by Predict *~ 1 9 and the actual solar irradiance in 
the k- th frame of the (d — l)-st day, denoted by G* -1 . Please 
note that we must also capture and predict the seasonal change 
of the sunny day solar irradiance, while filtering out random 
solar irradiance decay effects due to the presence of clouds. 
This is a smoothing operation. The Predict k value is 
calculated as follows 



P {Predict* 1 ,G* *) = 



Po> 



if Predict j?" 



1 < Gt 1 



P 0 e ^ Predlct k 1 G k 1 \ otherwise 

(18) 



where p o is the basis learning rate, and A is the decaying 
parameter for the learning rate. 

The motivation for this smoothing operation is as follows. 
Since (i) we want to predict the seasonal change of the sunny 
day solar irradiance while filtering out the effect of clouds, 
and (ii) G * _1 « Predict ^ _1 only if there are clouds, it is 
natural that our new predicted sunny day solar irradiance 
Predict * should not be so much influenced by G* ~ 1 (which is 
strongly affected by the clouds.) Therefore, we adopt the 
exponentially decaying learning rate (18), rather than the 
constant learning rate in the original exponential 
average-based prediction method. 

The intra-day refinement phase is performed at the start 
time of each k- th (for 1 < k < K ) frame in a day to predict 
the decay factors and thereby refining the prediction results 
over the rest of the frames in that day. The intra-day 
refinement is implemented as follows. Consider that we are 
currently at the start time of the k- th frame in the d-th day, i.e., 
T* k , and we intend to refine the prediction results of the solar 
irradiance in the k- th (for k < k < K) frame of that day. We 
denote the results of the intra-day refinement as Refine | (for 
k < k < K). Since at time T* k the actual solar irradiance in 
the (fc — l)-st frame is already known, we calculate Refine £ 
by 

Refine % «- (1 - y) ■ Predict % + y ■ ^ ■ Predict 

(19) 

where the coefficient 0 < y < 1. We learn the optimal y value 
by using the stochastic gradient descent method [25]. 
Typically, a y value in the range of [0.5,0. 7] will yield the 
best prediction results. 

The intuition for (19) is as follows: if the actual solar 
irradiance in the ( k — l)-st frame is higher than the predicted 
solar irradiance in that frame, i.e., G*_ t > Predict*_ ± , it is 
highly likely that the actual solar irradiance in the k -th 
(k < k < K) frame of the same day will also be higher than 
the predicted solar irradiance in that frame, and vice versa. 
Experimental results in Section VIII will demonstrate that the 
prediction error can be reduced to 42% of the initial prediction 
error by using intra-day refinement. 



Predict* = p{Predict* 1 ,G* x ) ■ G* 1 + (17) 

(1 - p (Predict* -1 , G * _1 )) ■ Predict*" 1 . 

In the above equation, the learning rate function 
P {Predict* -1 , G* -1 ) is set to 



VII. Global Control Algorithm 

We propose a cascaded feedback control based algorithm, 
comprised of an outer supervisory control loop and an inner 
control loop, to effectively solve the global control problem. 
The outer supervisor control loop is performed at each 
decision epoch (to be defined later) of the RTES-EH 
operational period, and it maintains the supercapacitor SoC 
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through effective task scheduling with DVFS, as well as 
selective task dropping. The inner control loop performs the 
optimal PV panel I-V operating point setting. The outer 
supervisory control loop and the inner control loop are 
performed every 300 seconds (which is equal to the frame 
length T'lcm ) and five milliseconds, respectively, in this paper. 

The decision epochs are defined as the start time of each 
frame. We denote the k- th (for 1 < k < K ) decision epoch by 
T sk , where subscript s stands for “start” and T sk = T sunrise + 
(k — 1)T lcm . At each decision epoch T sk , the supervisory 
control algorithm finds a schedule and frequency allocation of 
the task instances in the k- th frame, represented by fs t (t) (for 
1 < i < M and t E [T s k , T s k+1 ] .) The supervisory control 
algorithm is aware of future power harvesting (through solar 
irradiance prediction) and power consumption in the 
RTES-EH to perform effective supercapacitor SoC control. 

A. Motivation 

We have the following observations on the RTES-EH: 

1) At the beginning (i.e., in the morning) and at the end (i.e., 
in the evening) of the RTES-EH operational period, the 
solar irradiance level is low; while the solar irradiance 
becomes abundant in the middle (i.e., at noon) of system 
operational period in general. 

2) At the beginning of the RTES-EH operational period, the 
supercapacitor has no or very little energy — this is due to 
the effect of self-discharge overnight. 

3) The supercapacitor terminal voltage varies significantly as 
a function of the supercapacitor SoC. And the charger and 
DC-DC converter efficiencies are strongly dependent on 
the supercapacitor terminal voltage. 

Based on the above observations, we have the following 
motivations for the proposed global control algorithm. 

1) There exists a desirable amount of supercapacitor energy, 
such that the task instance drop rate can be minimized if 
the amount of energy stored in the supercapacitor is close 
to such a desirable value, due to the converter efficiency 
variation. The desirable energy value varies according to 
the solar irradiance change during a day. 

2) Selectively dropping some task instances at the beginning 
of the RTES-EH operational period may eventually help 
increase the total number of completely executed task 
instances since the energy stored in the supercapacitor may 
increase (towards the desirable value) along with system 
operation. 

3) It is desirable to use up the energy in the supercapacitor by 
the time T sunset , since the supercapacitor energy will be 
eventually dissipated due to self-discharge at night. 

B. Outer Supervisory Control Loop 

The objective of the supervisory control is to maintain the 
amount of supercapacitor energy (or equivalently, its SoC) 
around or above a desirable value, through effective task 
scheduling with DVFS as well as selective task instance 
dropping. In order to effectively perform supervisory control 



at each decision epoch, we first identify two new subproblems 
and provide corresponding solutions: the converter-aware task 
scheduling (CA-TS) problem (without task instance dropping) 
and the supercapacitor energy control (SEC) problem. Then 
we provide the supervisory control algorithm based on the 
solutions to the two subproblems, the above-described 
motivations, as well as effective solar irradiance prediction. 
Fig. 6 illustrates the relationship between the control 
algorithms and the subproblems in this paper in order for 
better understanding. As can be seen from Fig. 6, the 
previously defined SC problem and CA-FS problem are 
fundamental to the solutions to the CA-TS problem and the 
SEC problem, which are basics steps for the supervisory 
control loop. 




-> contains > based on the solution of 



Fig. 6. The relationship between control algorithms and subproblems. 

1) Converter- Aw are Task Scheduling (CA-TS) Problem 

Suppose we are at the k -th decision epoch T s k . The 
supercapacitor energy and terminal voltage are given by 
E cap(Ts,k ) an d V C ap(T s>k ), respectively. We are given a set of 
task instances, denoted by 3. The absolute release time and 
absolute deadline of all task instances in 3 are within the time 
period [T s k , T s k+1 ] (the k - th frame.) We are going to find a 
schedule and frequency allocation of all the task instances in 3, 
such that the amount of energy extracted from the 
supercapacitor is minimized. Please note that task instance 
dropping is not allowed in this problem. 

The method in [6] finds the optimal schedule and frequency 
allocation of a set of task instances on a continuously variable 
frequency processor, satisfying deadline constrains but 
without taking into account the converter power loss. On the 
other hand, the CA-FS problem proposed in Section V.B finds 
the optimal sensor node frequency for executing a specific 
task instance, taking into account the converter power loss but 
without any deadline constraint. Combining the method in [6] 
and the CA-FS problem, we can generate the optimal schedule 
and frequency allocation of a set of task instances satisfying 
deadline constraints and minimizing the energy extracted from 
the supercapacitor by considering the converter power loss. 
However, such allocated frequency levels are continuous 
variables. Then we employ the idea in [28] to use at most two 
discrete frequency levels for each task instance, and the two 
frequency levels are immediate neighbors of the continuous 
frequency level allocated previously. 
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Details of the proposed solution to the CA-TS problem are 
given in Algorithm 1. Similar to the solution to the CA-FS 
problem, we make the assumption in Algorithm 1 that the 
supercapacitor voltage will not change significantly during the 
time period [T s>k ,T sk+1 ] (the k - th frame). Based on this 
assumption, Algorithm 1 will provide the optimal discretely 
variable frequency allocation and schedule for all task 
instances considering the DC-DC converter power loss. 



Algorithm 1: Solution to the CA-TS Problem. 

Generate an optimal continuously variable frequency 
allocation and schedule for all task instances in 3, using the 
method in [6]. 

Find the optimal sensor node execution frequencies for task 
instances of T t , i.e., f° vt (Y ca p(T s k )) for 1 < i < M , by 
optimally solving the CA-FS problem. 

For each task instance: 

If the allocated frequency is less than f° vt (V cap (T s k )): 
Allocate the frequency f° pt (y c ap(T s>k )) to that task 
instance. 

Use the algorithm proposed in [28] to generate the optimal 
discretely variable frequency allocation and schedule for all 
task instances in 3. 



Next we calculate the estimated output energy of the 
supercapacitor. We derive the schedule and frequency 
allocation for task instances in the k- th frame, represented by 
fsft) for 1 < i < M and t E [T sk ,T sk+1 ] , by solving the 
CA-TS problem, in which we assume that the supercapacitor 
terminal voltage is V cap . We calculate the (estimated) 
supercapacitor discharging current / cap , O ut(0 based on 
fs(t) = Hi<i<M/ 5 i(0 and f Cflp , in the way described in 
Section IV. The estimated output energy of the supercapacitor 
during the k -th frame as a function of V cap , denoted by 
Eout.kiYcap), is calculated using 

Eout,k(Vca P ) = V cap ■ Sf +1 WO dL ( 21 ) 

Since the supercapacitor voltage is assumed to be a constant 
value V cap , the estimated output energy of the supercapacitor 
over each k -th frame ( k < k < K ) is also equal to 
Eout.k ( V C ap )• 

The estimated self-discharge energy loss in the 
supercapacitor from the /c-th frame to the /C -th frame as a 
function ofV cap , denoted by E sdktotal (y cap ), is calculated 
using: 



Return the schedule and frequency allocation. 



2) Supercapacitor Energy Control (SEC) Problem 
Suppose we are at the /c-th decision epoch T s k . The SEC 
problem aims to find the optimal amount of supercapacitor 
energy in the /c-th frame, denoted by E^ pk , such that the 
RTES-EH system availability will be maximized (i.e., the task 
instance drop rate will be minimized) when the amount of 
supercapacitor energy is near such a desirable value. In the 
following we discuss about the equivalent problem of finding 
the optimal supercapacitor terminal voltage V^ p k . 

For each possible supercapacitor terminal voltage value 
V cap , we calculate the estimated system availability as follows. 
For simplicity in estimation and without significant loss of 
accuracy, we assume that the supercapacitor terminal voltage 
will not change significantly (i.e., V cap (t) » V cap .) The 
estimated input energy of the supercapacitor from the /c-th 
frame to the K-th frame as a function ofl^ ap , denoted by 

Ein,k,totai(Vcap )> is calculated using 
E in,k,totaliYcap ) — 

Vcap ^T sk ^^^(^pi?(t)^pi?(t))Tcap(t) = l / cap ^cap,in(f') tit, (20) 

where the (estimated) supercapacitor input current / cap ^ n (t) is 
maximized by finding the optimal PV panel operating point 
(Vp V (t), I pv (t)) (i.e., optimally solving the SC problem) based 
on the predicted solar irradiance profile. Note that this is the 
software/model-based solution to the SC problem as described 
in Section V.A. 



^sd.k, total \Ycap) ^cap ' ^ cap, sd\Ycap) ' \J sunset i s,k ) ( 22 ) 

= (K-k + l)-C cap -^f-T LCM , 

where Icap.sdiVcap ) i s the supercapacitor self-discharge 
current as a function ofF cap . The above equation is derived 
from (8). 

We define the ratio of E in k total (Vcap) Esd.k, totally cap) t° 

E out.k (Vcap) as the estimated system availability. The 
estimated system availability represents the estimated number 
of frames that can be executed in the sensor node during the 
time period [T s k , T sunset ] . The optimal (target) supercapacitor 
terminal voltage Vcap,k the optimal V cap value that 
maximizes the estimated system availability. Therefore, Vla p k 
is calculated by 



Vcap.k = argmax^ 



£ in,k,total(ycap) @ sd,k,total(ycap) 
cap £out,k(V C ap) 



(23) 



3) Supervisory Control Algorithm 

The objective of the supervisory control algorithm is to 
maintain the amount of supercapacitor energy (or equivalently, 
its SoC) around or above a desirable value, through effective 
task scheduling with DVFS as well as selective task instance 
dropping. Consider the supervisory control algorithm 
performed at the decision epoch T s k . The supercapacitor 
energy and terminal voltage are given by E cap (T sk ) and 
Vap(Ts,k )> respectively. The optimal (target) supercapacitor 
energy is given by E% r pk , satisfying E^ pk = \c cap (V ^ ik ) 2 . 
We use e cap k = E cap (T s k ) — E^ p k to denote the difference 
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between the amount of supercapacitor energy and the target 
value. The supervisory control algorithm goes as follows: 

Case I (e capk > 0): In this case the amount of energy stored 
in the supercapacitor is abundant (higher than the target value.) 
We schedule all the task instances in the k- th frame using 
Algorithm 1, without task instance dropping. 

Case II ( e cap,k < 0): In this case the objective of the 
supervisory control algorithm is to make the supercapacitor 
energy at the next decision epoch, E cap (T sk+1 ), equal to or 
higher than the value E cap (T sk ) + 5 ■ \e capk \, in which 6 is a 
predefined feedback control parameter. We may have to drop 
some task instances to achieve this goal. The basic idea is to 
keep dropping the most energy consuming task instance and 
performing task rescheduling using Algorithm 1, until the 
estimated amount of energy stored in the supercapacitor at 
time T sk+1 satisfies the requirement. Detailed procedure is 
shown in Algorithm 2. 

However, there is no need for keeping the supercapacitor 
energy around or above the target value near the end of the 
RTES-EH operational period (i.e., in the evening.) We can 
schedule all the task instances without dropping in this case. 
We use the following estimation-based procedure to decide at 
a decision epoch T s k whether performing supercapacitor 
energy maintenance is necessary. We first calculate the 
estimated supercapacitor input energy E inXtotal (V cap (T sk )) 
and self-discharge energy loss E sdXtotal (y cap (T sk )) using the 
method described in Section VII.B.2). Next we calculate the 
estimated supercapacitor output energy E outk total (V cap (T sk )) 
as 

^out,k, total ( Vcap(js,k )) = (K ~ k + 1 )E outk (V cap (T sk )), (24) 
where E out k (v ca p(js,k)^) is calculated in the way described in 

( 21 )* If ^in,k, total (^cap CTs,/c)) apis's, k) — 

E out, k, total (fcap(^5,/c)) ^sd,k, total (V cap (T Stk )), we conclude 

that the supercapacitor energy plus the future harvested energy 
is enough for scheduling all the task instances until T sunset , 
and hence the supervisory control algorithm will schedule all 
the task instances in the k- th frame without task instance 
dropping. Otherwise, we perform supercapacitor energy 
maintenance based on selective task instance dropping. Details 
of the proposed supervisory control algorithm are given in 
Algorithm 2. We again make the assumption that the 
supercapacitor terminal voltage will not change significantly 
during the k-th frame. 



Algorithm 2: The Supervisory Control Algorithm. 

Predict the solar irradiance levels of the rest of the frames in 
the day, i.e., if k= 1, calculate Predict f for k < k < K ; 
otherwise, calculate Refine | for k < k < K . 

Calculate the Ei n k)tota i (Vcapi^s.k)^ E out, k, total (lcap(Ts,/c))> 



and E sd)ktota i ^J^ a pCTs,/e))- 

Initialize the set 3 to include all task instances in the k - th 
frame. 

If Ein,k, total capias, k )) E capias, k ) — 

E out, k, total (Vcapfts,k)) ^sd,k, total 

Use Algorithm 1 to generate a schedule for all task 
instances in 3. 

Else: 

Calculate E^ p k by solving the SEC problem. 

Calculate e caPik = E cap (T sk ) - E^ k . 

If e cap ,k ^ 0 : 

Use Algorithm 1 to generate a schedule of all task 
instances in 3. 

Else: 

Calculate the estimated supercapacitor input energy 
during the k-th frame, given by Ein,k ( V C ap{T s ,k )) — 
Vcap(^s,k) ' 

^T sk ^^^(V r p V (t) ; /p V (t));V r ca p(t)=V r ca p(7 7 s k) I cap, in (0 tit. 

Calculate the estimated supercapacitor self-discharge 
energy loss during the k-th frame, given by 

p (v (T W - C ( WCU,fc )) 

C'sd.k \fcap \ l s,k)j L cap r l LCM- 

Use Algorithm 1 to generate a schedule of all task 
instances in 3. 

Calculate the estimated supercapacitor output energy 
E ou t,k ( Vcap(js,k )) during the £-th frame based on the 
schedule. 

E capias, k) Et nk (V C ap(fs,k )) — 

Eout,k (Vcap(Ts,k )) — ^ sd,k (Vcapfts.k)^ ^ ^capiTs.k) 

"f ^ ' I t?cap,k I • 

Drop the most energy consuming task instance from 

3. 

Use Algorithm 1 to generate a schedule for the new 
set 3. 

Calculate E outk (Vcap(Ts,k)) based on the new 
schedule. 

Return the schedule and frequency allocation. 



C. Inner Control Loop 

The inner control loop determines the optimal I-V operating 
point of the PV panel (and hence the input current of the 
supercapacitor). The inner control loop is performed every 
five milliseconds in this paper. At time t, the supercapacitor 
energy is given by E cap {f). The inner control loop determines 
the optimal PV panel operating point (V pv (f),I pv (t)) by 
optimally solving the SC problem based on the actual solar 
irradiance C(t). Please note that this is the hardware-based 
solution to the SC problem as described in Section V.A. 
Besides, the sensor node execution frequency fs(t) is given 
by the task schedule and frequency allocation generated by the 
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supervisory control algorithm. The inner control loop also 
controls the sensor node supply voltage V proc (t) using the 
DC-DC converter, so that the embedded sensor node properly 
runs at frequency fs(t ). 

VIII. Experimental Results 

First, we need to test the accuracy of our solar irradiance 
prediction algorithm. We use the solar irradiance profiles 
measured at Duffield, VA, for the year of 2007. Remember 
that the initial prediction for the sunny day solar irradiance is 
based on the actual solar irradiance in the previous day, and 
therefore, we use the solar irradiance profiles for the first 91 
days in the year as the training data and we predict the solar 
irradiance for the rest of days in the year. 

Fig. 7 shows the actual and predicted solar irradiance 
profiles of day 170, 230, 290, and 350 at Duffield, VA in the 
year of 2007. Please note that the solar irradiance levels are 
normalized by the solar irradiance value G STC at standard test 
condition. Fig. 7 contains the actual solar irradiance profiles, 
the initial prediction result, and the refined prediction result. 
Remember that the initial prediction is performed at the 
beginning of each day (i.e., 6:00am in this work) and the 
intra-day refinement is performed at the start time of each k - th 
(for 1 < k < K ) frame in a day. In Fig. 7, the refined 
prediction result for the k - th frame in a day is obtained (and 
refined for the last time) at the start time of the k - th frame. We 
can observe from Fig. 7 that (i) in a sunny day (e.g., day 290) 
the initial prediction is accurate enough and the intra-day 
refinement cannot improve accuracy further and (ii) in a 
“cloudy” day the intra-day refinement improves prediction 
accuracy significantly. On average, the prediction error of 
initial prediction is 19.4% and the prediction error of refined 
prediction is 8.1%. 




Fig. 7. The actual and predicted solar irradiance profiles of four different days 
(day 170, 230, 290, and 350) at Duffield, VA in the year of 2007. 



In the next, we compare our proposed global control 
algorithm with a baseline control algorithm on the RTES-EH 



in terms of task instance drop rate. The PV panel used in the 
RTES-EH has a maximum output current of 0.305 A, a 
maximum output voltage of 27.8 V, and an MPP output power 
of 6.2 W under the standard test condition. The I-V and P-V 
characteristics of the PV panel are shown in Fig. 5. We use the 
Finear Technology FTM4607 converter as the charger and 
voltage converter model in the RTES-EH. A 500 F 
supercapacitor is used as the energy storage device. The 
embedded sensor node in the RTES-EH has six discrete 
operating frequencies: 500 MHz, 600 MHz, 700 MHz, 800 
MHz, 900 MHz, and 1000 MHz. Correspondingly, the sensor 
node has six discrete supply voltage levels: 1.0 V, 1.2 V, 1.4 

V, 1.6 V, 1.8 V, and 2.0 V, and six power levels: 0.62 W, 0.89 

W, 1.21 W, 1.58 W, 2.00 W, and 2.47 W. The task set 
executed in the RTES-EH consists of five tasks i.e., 
T lt T 2 r" the periods D t (1 < i < 5) of which are 50 s, 60 
s, 100 s, 150 s, and 300 s, respectively. The frame length T LC m, 
which is the hyper-period of all the tasks, is therefore 300 s in 
the RTES-EH. The task instance execution times i.e., 
WJ f r eq ma x (1 ^ t < 5), which is the execution time of an 
task instance of T t when the sensor node runs at its maximum 
operating frequency freq max , are 10 s, 10 s, 15 s, 15 s, and 30 
s, respectively. The above task setup is based on habitat 
monitoring sensor node [26] and wireless camera sensor node 

[27] . The baseline control algorithm is based on the same 
RTES-EH system architecture. The baseline algorithm adopts 

[28] to find a discretely variable frequency allocation and 
schedule of all task instances in the sensor node, and 
incorporates a simple procedure to protect the supercapacitor 
from energy depletion. When the sensor node is about to 
execute a new task instance, the baseline controller checks 
whether the current supercapacitor terminal voltage is less 
than a predefined threshold value. If so, that task instance will 
be dropped. 




Fig. 8. Task instance drop rate from the proposed global control algorithm and 
baseline control algorithm based on solar irradiance profile of day 170. 

Fig. 8 and Fig. 9 demonstrate the task instance drop rate 
from the proposed global control algorithm and the baseline 
control algorithm when the solar irradiance profiles of day 170 
and day 290 are used, respectively. The task instance drop rate 
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for different supercapacitor initial voltage V ca p (Tsunrise ) levels 
(i.e., from 1 V to 4 V) are plotted in Fig. 8 and Fig. 9. Please 
note that the proposed global control algorithm is based on the 
solar irradiance prediction results, whereas the baseline 
control algorithm is based on instantaneous solar irradiance 
level. We can observe that the proposed global control 
algorithm consistently outperforms the baseline control 
algorithm. The proposed global control algorithm achieves a 
maximum of 63% reduction in task instance drop rate 
compared with the baseline algorithm. 




Fig. 9. Task instance drop rate from the proposed global control algorithm and 
baseline control algorithm based on solar irradiance profile of day 290. 




Fig. 10. Supercapacitor terminal voltage change from T sunrise to T sunset (i.e., 
a system operational period). 



Furthermore, we examine the supercapacitor terminal 
voltage V cap (t ) change and number of dropped task instances 
during a system operational period (from T sunrise to T sunset ) 
shown in Fig. 10 and Fig. 11, respectively. Combining the 
information in Fig. 10 and Fig. 11, we can find that at the 
beginning of the system operational period the baseline 
control algorithm tends to execute all the task instances until 
V cap (t) drops to a threshold value, whereas the proposed 
global control algorithm is more conservative by dropping 
some task instances even though the supercapacitor has 
enough initial energy. After some time, the supercapacitor in 
the global control algorithm can accumulate higher energy 
than the supercapacitor in the baseline control algorithm when 
the solar irradiance is abundant (at noon). Therefore, the 
global control algorithm can maintain a lower task instance 



drop rate in the late afternoon (when the solar irradiance is not 
abundant). 




Fig. 1 1 . Number of dropped task instances per frame length Tlcm from 
Tsunrise to T sunset (i.e., a system operational period). 

In addition, we demonstrate the difference between the task 
instance drop rates from the proposed global control algorithm 
when using the predicted and actual solar irradiance profile in 
Fig. 12. The solar irradiance profile of day 230 in Fig. 7 is 
used for the simulation. There are two observations from Fig. 
12. First, the general task instance drop rates are higher 
compared with the results in Fig. 8 and Fig. 9. This is due to 
the relatively low solar irradiance in day 230. Second, the task 
instance drop rate from the proposed global control algorithm 
using predicted solar irradiance profile is slightly higher than 
that using actual solar irradiance profile. It demonstrates that 
our proposed global control algorithm is robust under solar 
irradiance prediction error. 




Supercapacitor Initial Voltage (V) 



Fig. 12. Task instance drop rate from the proposed global control algorithm 
using predicted and actual solar irradiance profile of day 230. 



IX. Conclusion 

In this paper, we investigate the question of how one can 
achieve the prolonged lifespan of a real-time embedded 
system with energy harvesting capability by applying a global 
control strategy. The global controller performs optimal 
operating point tracking for the PV panel, SoC management 
for the supercapacitor, and energy-harvesting- aware real-time 
task scheduling with DVFS for the sensor node, based on an 
accurate solar irradiance prediction method. The controller, 
which accounts for dynamic I-V characteristics of the PV 
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panel, terminal voltage variation and self-discharge of the 
supercapacitor, and power losses in power converters, 
employs a cascaded feedback control structure with an inner 
control loop determining the I-V operating point of the PV 
panel and an outer supervisory control loop performing 
real-time task scheduling and setting the voltage and 
frequency levels in the sensor node. 
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